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Gravitational-wave searches for the merger of compact binaries use matched-filtering as the 
method of detecting signals and estimating parameters. Such searches construct a fine mesh of 
filters covering a signal parameter space at high density. Previously it has been shown that singular 
value decomposition can reduce the efi'ective number of filters required to search the data. Here we 
study how the basis provided by the singular value decomposition changes dimension as a function 
of template bank density. We will demonstrate that it is sufficient to use the basis provided by 
the singular value decomposition of a low density bank to accurately reconstruct arbitrary points 
within the boundaries of the template bank. Since this technique is purely numerical it may have 
applications to interpolating the space of numerical relativity waveforms. 



I. INTRODUCTION 

Several broadband laser interferometer gravitational- 
wave (GW) detectors are operating at high sensitivities 
and will continue to improve over the next decade [ - ] . 
As detectors improve it is increasingly likely that GW 
astronomers will observe gravitational radiation emitted 
from the coalescence of compact binary systems involving 
neutron stars and or stellar mass black holes [' ] . 

Because compact binary coalescence (CBC) waveforms 
are well modeled, GW searches for such signals are con- 
ducted by matched filtering the detectors' data with 
banks of template waveforms, chosen to adequately cover 
a region of the signal parameter space [ ]. For GW sig- 
nals from the merger of compact objects with negligi- 
ble spin, this parameter space is defined by functions 
of the masses of the two objects. To search for sig- 
nals within this parameter space, a bank of templates 
is constructed to sample the parameter space sufficiently 
densely such that there is minimal loss of signal-to- 
noise ratio (SNR). Traditionally, template banks used 
to search this two-dimensional signal parameter space 
have been constructed using the (^2)* lattice [8], referred 
to as "hexagonally-placed" template banks. This prob- 
lem becomes more difficult in higher dimensions, where 
other types of template placement algorithms have re- 
cently been investigated [-*-! "-]. 

In [ ] the singular value decomposition (SVD) was ap- 
plied to CBC waveforms to show how hexagonally-placed 
template banks with M templates could be implemented 
with N' <^ 2M filters (2M being the nominal number of 
filters required for the M 2-phase templates). This was 
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achieved by truncating the SVD of the matrix consisting 
of the time-series of the template waveforms. Here we 
demonstrate that the bases identified by the SVD is effec- 
tive at spanning the space of all CBC waveforms within 
the region of parameter space sampled by the original 
bank. We find that the SVD of a low-density bank pro- 
vides a basis suitable for constructing all the waveforms 
from a higher-density bank, even waveforms at arbitrary 
locations within that region of parameter space. 

This paper is organized as follows. Sec. II describes 
how we apply the SVD to approximately embed the sig- 
nal manifold in a vector space. Sec. Ill tests this em- 
bedding by reconstructing various points in the mani- 
fold. Finally, Sec. IV discusses possible applications of 
this technique. 

II. ENCLOSING THE SIGNAL SPACE WITH 
SINGULAR VALUE DECOMPOSITION 

In this section we explore how the number of basis 
vectors required to reconstruct a template bank scales 
with the initial density of the template bank. We define 
a template bank of signal waveforms covering a patch 
V of the signal manifold, which is used to test for the 
presence and strength of signals from V in the detec- 
tors' data. We construct a signal matrix in the same 
manner as [13]. Specifically, we create a real- valued ma- 
trix H by alternately filling its rows with the real and 
imaginary parts (cosine and sine) of the template wave- 
form time series from a CBC template bank covering V, 

H = = {mi,^hum2,^h2, ...MmMmV ■ 

As in [ ], we constructed the template matrix with 
chirp masses = Mr]^^^, where M = 7711+7712 is the to- 
tal mass and rj = mim2/M^ is the symmetric mass ratio, 
of 1.125Af0 < Mc < I.24OM0 and component masses of 
IMq < 7711,7(72 < 3Mq. Template banks covering this 
region are created using template placement algorithms 
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of the LIGO Algorithms Library [14]. Template place- 
ment is done in the (to,T3) plane, where tq and T3 are 
defined as 

ro^i-,infor'^'M-'/\ (1) 
r3 = f (^/o)''/'M-%-3/5, (2) 

and where /o is some fiducial frequency, which we choose 
to be /o = 60 Hz. 

The non-spinning waveforms for each template are pro- 
duced to 3.5 post-Newtonian (PN) order, sampled at 
2048 Hz, up to the Nyquist frequency of 1024 Hz. The 
last 10 seconds of each waveform, whitened with the ini- 
tial LIGO amplitude spectral density, are used to con- 
struct H. The SVD is then applied to H, decomposing 
the matrix into two unitary matrices, V and U, and a 
diagonal matrix S 

H = VSU"^, (3) 

where U is a matrix composed of basis vectors (i.e., unit- 
norm time-series vectors), V is a matrix composed of 
reconstruction coefficients, and S is a matrix containing 
the singular values of H. 

In [13], it was demonstrated that truncating the re- 
construction of H to use only the N' basis vectors with 
the largest singular values results in an average fractional 
SNR loss (Sp/p) proportional to the sum of the discarded 
singular values squared. In this investigation, we trun- 
cate these reconstruction matrices at {Sp/p) = 10"''. 
This corresponds roughly to the truncation error of IEEE 
754 32-bit floating-point numbers. 

We explore how the number of basis vectors changes as 
the number of rows in H is increased by generating tem- 
plate banks for V with increasing density (i.e., increasing 
minimal match) . We confirmed that the number of basis 
vectors required to reconstruct H saturates at a partic- 
ular value of minimal match. Fig. 1 shows that as the 
minimal match of the template bank is increased, result- 
ing in denser samplings of V, the number of basis vectors 
needed to reconstruct H to the required accuracy satu- 
rates around a minimal match of ~ 89.9%. This indicates 
that V is able to be embedded — to an accuracy of 1 part 
in 10^ — in a vector space consisting of ~ 150 dimensions. 

In the next section we will demonstrate how the basis 
waveforms identified by the coarsely sampled bank can 
be used to reconstruct templates at arbitrary points on 
the signal manifold. 

III. EFFICIENT RECONSTRUCTION OF 
WAVEFORMS IN THE MANIFOLD 

In order to determine how well these waveforms can 
be reconstructed, we compute a quantity called the av- 
erage fractional SNR loss 5pa/pa- This quantity can be 
thought of as the mismatch between the original wave- 
form ha and the reconstructed waveform /i^, averaged 
over the phase angle. It tells us how far the reconstructed 
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FIG. 1. The number of filters as a function of minimal match, 
which increases with the density of the template bank. The 
total number of filters in the template bank. A'', is shown by 
the dashed line. The number of filters needed to reconstruct 
the template matrix such that {5p/ p) — , N' , is shown 
by the solid line. We find that the number filters needed 
to reconstruct H saturates when the minimal match reaches 
~ 89.9%. 




FIG. 2. A visual representation showing the two types of 
points of V that we can choose to reconstruct. The left car- 
toon shows an example point that is in H, and thus in V. 
The right cartoon shows an example point that is not part of 
H but is in V. 

waveform is from the original waveform. As in Eq. (25) 
of [ ], (5/9q//9q, is given in terms of SVD quantities as 

Sp 1 ^ 

where W(2q-i)^ and W(2a)/i are the reconstruction coeffi- 
cients for the real and imaginary parts, respectively, of 
the ath waveform associated with the pth basis vector 
and are elements of V, cr^ is the pth element of S, and 
the sum is over the truncated terms of V and S. 

We test this embedding of the signal manifold to see 
how well various points in the manifold can be recon- 
structed. The tests points we reconstruct are of two 
types: 1) those from the original signal matrix H, 2) 
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are then reconstructed using 



h' 




(6) 
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FIG. 3. A histogram of the reconstruction accuracy of the 
waveforms that went into the construction of H. As expected, 
the average reconstruction (solid hne) accuracy matches the 
expected fractional SNR loss (dashed line). 



those absent from H but within V. These two types of 
tests are illustrated in Fig. 2. 

A test of the first type is shown in Fig. 3. This shows 
that the average reconstruction accuracy for points from 
H agrees with our chosen value of 10"''. This result is 
expected as it is an extension of the investigation from 
Fig. 4 of [13] applied to a more stringent reconstruction 
accuracy. 

A test of the second type is shown in Fig. 4. To choose 
points uniformly from V but absent from H, we generate 
a denser template bank within the same region of param- 
eter space described in the Sec. II. Specifically, we gener- 
ate this template bank with a minimal match of 99%. In 
order to test the reconstruction accuracy of these wave- 
forms, we project the real and imaginary parts of the 
waveforms onto the basis vectors from the SVD of H 
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(5) 



where v'^^ represents a reconstruction coefficient associ- 
ated with the /ith basis vector for the real or imaginary 
part of the ath waveform from the denser template bank, 
(Tp is the /ith element of S, haj is the jth time sample 
of the real or imaginary part of the ath waveform from 
the denser template bank, and is the jth time sample 
from the /ith basis vector in U. The real and imaginary 
parts of the waveforms from the denser template bank 



The distribution of Spa/pa for these waveforms, left 
panel of Fig. 4, shows a tail extending to large mis- 
matches. Examining where these large mismatches are 
located in parameter space, we find they originate from 
near the boundaries of V. Removing the test points near 
the boundaries in the tq direction, shown in Fig. 5, we 
find the tail of large mismatches disappears. 

An additional test of the second type, which systemat- 
ically explores the reconstruction accuracy near a point 
whose waveform went into H, is shown in Fig. 6. The left 
panel shows a set of three nearest-neighbor templates. 
We investigate how the reconstruction accuracy varies 
as one moves from point A to the central point, point 
B. Point B is assumed to have the largest mismatch 
between its waveform and the waveforms from any of 
the three surrounding points. We also compute the mis- 
match between the waveforms along AB and the wave- 
form of point A with and without maximizing over phase 
and time, called the fitting factor and match respectively. 
The fitting factor falls to the minimal match of the tem- 
plate bank when comparing the waveforms from A and B, 
which is expected as the minimal match involves maxi- 
mizing over phase and time. The reconstruction accuracy 
associated with SVD projection is consistently high and 
close to the chosen reconstruction accuracy of 1 part in 
10^. 



IV. DISCUSSION 

These investigations show that the SVD can be used 
to find a set of basis vectors that not only span the signal 
matrix H, but also enclose the signal manifold V sampled 
by H. 

GW pipelines that search for known waveforms, such 
as GWs from CBCs, commonly compute waveform con- 
sistency statistics that compare the observed response of 
a template waveform filter to the data with what one 
would expect given the presence of that signal. These 
consistency statistics are found to perform better when 
the mismatch between the template waveform and the 
signal waveform is small [ ' "]. Filtering with a fixed den- 
sity template bank can introduce mismatch between the 
nearest template and the signal. This mismatch can be 
greatly reduced if one is able to find the exact point in 
parameter space where the signal is located and filter the 
data using that point. Using the SVD basis vectors, one 
could reconstruct a point closer to the point of the signal 
and improve the waveform consistency statistics. 

Parameter estimation techniques for GWs from CBCs 
often use Monte Carlo Markov Chain algorithms to 
search the parameter space. This involves producing 
waveforms and filtering the data against many points of 
the parameter space. The SVD could also be used to 
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FIG. 4. {left) A histogram of the mismatch between waveforms from V that were not in H and their reconstructed versions. 
The peak of the mismatch is slightly above the expected fractional SNR loss for waveforms from H (dashed line), {right) How 
the mismatch varies across V, averaged across the tz direction. The largest mismatches come from near the borders of the 
template bank in the tq direction. Fig. 5 restricts our attention to the central 75% of the domain of V, whose boundaries are 
shown as vertical dotted lines. 
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FIG. 5. Similar plots as in Fig. 4, eliminating test points 
near the boundaries in the to direction, (left) A histogram 
of the mismatch between waveforms from V that were not in 
H and their reconstructed versions. The average reconstruc- 
tion (solid line) accuracy is slightly worse than the expected 
fractional SNR loss (dashed line). 



interpolate waveforms that are expensive to compute, as 
in the case of waveforms produced by solving differential 
equations. Also, if one filtered the data using the basis 
vectors from the SVD, it would be very easy to recon- 
struct to high accuracy the output one would see if one 
had filtered the data using a waveform from anywhere 
within the parameter space. 

In order to gain benefit from these applications, it 
would be necessary to determine the reconstruction coef- 
ficients in a computationally efficient manner. This work 
has not tried to address this problem as: 1) it has as- 
sumed the target waveforms are known, and 2) it com- 
putes the reconstruction coefficients using computation- 
ally expensive inner products. Generation of these recon- 
struction coefficients warrants future investigation as the 
benefits derived from this technique would be substan- 
tial. 
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FIG. 6. (left) A plot showing a smaller region of V. The circles are points whose waveforms go into H. The line segment AB 
connects one of those points, A, with the point B, the central point of A and two of its nearest neighbors. Point B is situated 
such that it should have a 89.9% fitting factor with each of the surrounding points, [right) How well a waveform from a given 
point of AB, hp, can be "matched". The Match shows the normalized inner-product between /ia and the waveform from the 
corresponding point along AB. The Fitting Factor shows the same match maximized over phase and time. Since we are using 
a 89.9% minimal match bank, it is expected that the fitting factor falls to around that value. The SVD Projection is given by 
1 — 5 Pa/ pa- The SVD basis vectors are able to reconstruct to high accuracy all points along the line. 
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